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In this work we consider a problem related to the equilibrium 
statistical mechanics of spin glasses, namely the study of the Gibbs 
measure of the random energy model. For solving this problem, new 
results of independent interest on sums of spacings for i.i.d. Gaussian 
["V 1 ^ random variables are presented. Then we give a precise description 

p I ■ of the support of the Gibbs measure below the critical temperature. 

, 1. Introduction. The problem of spin glasses is considered a challenge 

for the 21st century from both the theoretical physics point of view [47] 
and the mathematical point of view [46]. Spin glasses are alloys like Au-Fe 
that have a very small density of Fe. At very low temperature they present 
^t. • remarkable magnetic properties. On the experimental level, measurements 

of the magnetic susceptibility of such alloys at very low temperature were 
done by Cannela and Mydosh [8]. 

As mentioned in [8] , the first model used to describe such an alloy was the 
, Rudderman, Kittel, Kasuya and Yosida (RKKY) model, a long-range model 

^£ ' with alternating sign in the interaction between the Fe atoms. Keeping just 

the alternating sign of the interaction, Edwards and Anderson introduced 
' ■^ | a model with short-range random coupling [21]. Despite a lack of rigor- 

ous results for the Edwards-Anderson model, a lot of theoretical progress 
has been made (see [36, 37, 48] and references therein). As a caricature 
of models with a spatial structure like the Ising model, mean field mod- 
els such as the Curie-Weiss model have been introduced to give a simple 
partial description of physical phenomena, namely the existence of sponta- 
neous magnetization at low temperature. The mean field model associated to 
the Edwards-Anderson model is called the Sherrington-Kirkpatrick (S-K) 
model [45]. From the theoretical physics point of view, much work had been 
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completed on the S-K model before Parisi [39, 40] introduced his famous 
replica symmetry breaking argument. A very important work by Mezard, 
Parisi and Virasoro [34] followed (see also [35]). Rigorous results on the S-K 
model can be found in [[1, 2, 10]-[12, 25, 28, 31, 46]]; however, we are still 
far from having complete results and this subject remains a very active field 
of research. 

Due to the difficulty of studying the S-K model, even from a nonrigorous 
point of view, various mean field models have been introduced. Among them 
is the random energy model (REM) [16, 17], which is considered to be the 
simplest model of spin glasses (see [27]), and a whole class of models called 
generalized random energy models (GREM) that present a tree structure 
[[18]-[20]]. Basic rigorous results on the existence of the free energy for the 
REM and its fluctuations are given in [7, 22, 26, 38], whereas for the GREM, 
see [[9], [26]] and, more recently, [[5, 6]]. A very important fact is that the 
free energy of the GREM can also be found by using the replica symmetry 
breaking argument of Parisi (see [20] ) , since we know rigorously the explicit 
expression of the free energy (see [9]); this is an important test for the 
validity of the Parisi theory. Various pedagogical expositions on the REM 
and the GREM are available (e.g., [3, 4, 42, 44]). 

The Hamiltonian or energy function of the REM can be written as 

ATl/2 

(1-1) #REM(0") = ~1^H 2-j Ja&a, 

aC{l,...,N} 

where the sum is over all the 2^ subsets a of {1, ...,N}, (J a ,a C {1, . . . ,N}) 
is a family of i.i.d. standard Gaussians defined on a common probability 
space (CI, £,P) and a a = JJ iea o% with a = 1. The above way of writing the 
Hamiltonian of the REM emphasizes the mean field aspect of the model, 
since all spins interact with a strength that does not depend on their dis- 
tances. 

It is easy to see that the Hamiltonians (Frem(c), cr £ { — 1 5 +1} ) of the 
REM form a family of 2 N i.i.d. Gaussian random variables with mean and 
variance N, defined on (fi,£,P); see [24]. In particular, we can also write 
(H RE M(a),ae {-1, +1}") = (-VNX a ,a€ {-1,+1}*), where (X a , a G {-1, +1}*) 
is a family of i.i.d. Gaussian random variables with mean and variance 1, 
or even (—y/NXi, 1 < i < 2 N ) when the relationship with the spins a is ir- 
relevant. Here y/N ensures that the energy in a volume N is P-almost surely 
of order N, as required by statistical mechanics theory. 

Given (3 > 0, the inverse temperature, let us denote by 



(1.2) 



Z N = Z N (I3)= J2 e- 1311 ** 1 "^ 

ae{-l, + l} N 
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the finite volume partition function and by 
(1-3) F N (J3) = -^log Z N (0) 

the finite volume free energy. 

It was proved in [38] that for all f3 > 0, limjv^oo F N ([3) = F(/3) exists P- 
almost surely and in L P (Q,¥) for 1 < p < oo (note that what is called free 
energy in [38] we denote —(3Fn (/?)). The nonrandom function F((3) is twice 
differentiable in 0, whose second derivative has a jump at (3 C := v / 2Tog~2. In 
fact, F((3) is equal to -(3/2 - /3 C 2 / '(2/3) for j3 < j3 c and -f3 c for (3 > C , as 
expected from the results of [16]. In the physics literature this is called a 
third order phase transition, as mentioned explicitly in the title of the first 
rigorous paper on the REM (see [22]). More results on the asymptotics of 
F N (J3) can be found in [7, 26, 38]. 

What is important to note is that when (3 > (3 C , the main contribution 
to Fn((3) comes from the terms in Z^{f3) that have the lowest possible en- 
ergies, which are the a that have -£7rem(°" ) of order — iV-y/21og2 = —N(3 C , 
that is, the (P-almost sure) asymptotic value of the minimum of 2 N Gaus- 
sian random variables with mean and variance N; this fact is explicitly 
mentioned in all pedagogical expositions (e.g., [3, 4, 42]). In the physics lit- 
erature, such a system is said to be frozen in the sense that, in the whole 
range of (3 > [3 C , the lowest possible energies give the main contribution to 
the free energy, while in a "nonfrozen" system, this contribution occurs only 
at zero temperature. This can be seen on F(f3) since when (3 > (3 C , the free 
energy does not depend on (3. 

To study the fluctuations of the free energy Fn {(3) , we consider the num- 
ber of random variables of the sample that are below some well chosen 
nonrandom energy, and then use classical convergence to a Poisson point 
process [30] to describe some fluctuations of the model. This was done in 
[26] for the Boltzmann factor exp(— (3Hrem(&)) (see also [3]) and in [7] for 
the partition function. 

The aim of this article is to study the finite volume Gibbs measure hn,/3 
of the REM when (3 > f3 c and N is large enough. Here Hn,p is defined for 
each ./V as the random probability measure on {-1,+!}^ which gives the 
configuration a the weight 

(i-4) mA a ) = — 7 > 

where is defined as in (1.2). 

For a given sample (^remW^ £ {— 1, +1}- ), when f3 > (3 C , the main 
question is: What are the sample: dependent configurations a, where the 
Gibbs measure is concentrated, when N is very large? 
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Note first that the sample-dependent configuration which corresponds 
to the minimal value of the sample -Hrem (f) is clearly among these configu- 
rations. So instead of considering the variables that are below a nonrandom 
energy as was done to get the Poisson point process mentioned above, it is 
better to consider the variables that are above the minimal one that is a 
sample-dependent energy. It is clear that the order statistics of the random 
variables -£Trem(c)> namely 

(1.5) H REM (a^) > H REM (^ 2N ^) > ■ ■■> #rem(<7 (2) ) > #remG7 (1) ), 

come into play. To obtain information on the support of the Gibbs mea- 
sure, we can subtract the minimal energy from all the Hamiltonians; there- 
fore, the spacings -ffREM(c^ fc+1 ^) — -Hrem^^) and the sums of spacings 
-f^REM^^)) — -ffREM(c^) are the basic objects we need to study. 

Since our main objective is more the study of the Gibbs measure of the 
REM than the sum of spacings per se, we first present a list of questions 
that we asked ourselves when we started this work. Then the strategy we 
used and the results stated in the next section about the sums of spacings 
will be clearer. 

First, how many successive terms k = k^ in the sum of spacings do we 
need to take to have a "good" approximation of the Gibbs measure when 
considering the probability measure on { — 1,-1-1}^ that havs only these k 
terms? 

There are two points to be answered precisely: What does "good ap- 
proximation" mean and in what P-probability sense can we expect that 
such approximation holds? For the first point we choose the total varia- 
tion distance for probability measures on {— Concerning the P- 
probabilistic sense, some care is needed. We first note that there is a scaling 
factor ^/N in the definition of the Hamiltonian, as mentioned previously. 
Recalling that we have a sample of 2 N Gaussian random variables, this 
factor can be written, up to a constant, as V2\ogW. So the first ques- 
tion is related to the behavior of V 2 log 2 N [//rem (c^ ) — -£7rem(c^)] as a 
function of k for k = k^ that diverges with N . Since we have chosen the 
total variation distance, it will be sufficient to neglect all the terms in the 
Gibbs measures that are larger than k^. An easy way to do this is to re- 
quire that Y%Lk N exp -{{(3 / (3 c ){^2\og2 N [H KEM {a^) - H REM (a^)})} goes 
to zero. The strongest P-probability sense we can expect to get is P-almost 
surely. At this point, it is important to note that we do not expect that 
the Gibbs measure itself, or its approximation converges almost surely since 
they merely converge in law. However, their total variation distance, being 
a difference, could perfectly converge almost surely to 0. A similar fact was 
proved not for the Gibbs measure as a measure on {—1,-1-1}^, but for an 
induced measure in the random field Curie- Weiss model [? ] . The worst that 
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could happen by taking such a strong convergence would be that we get a 
/cat that is rather large. 

The main result we need for the sum of spacings is given in Proposition 
2.2. It is stated for a sample of n Gaussian random variables, instead of 2^, 
since these results have independent interest. 

The second question is: Is the approximated Gibbs measure, that is, the 
one with only k^ terms obtained previously, the uniform measure on k^ 
points chosen without replacement within 2 N points, a point mass at the 
minimum or some other value? To be able to answer such questions, the 
relative weights of the approximated Gibbs measure and, therefore, the first 
k^ sum of spacings come into play. We use a classical representation of 
the Gaussian random variables in term of uniform random variables, mainly 
because a lot of explicit distributions for the spacings of uniform random 
variables are available. The results needed for the sum of spacings are stated 
in Proposition 2.3, where the first k^ sums of spacings can be represented 
as successive partial weighted sums of exponential random variables. 

Note also that since the kjy was found to be able to neglect some tails of 
the Gibbs measure, P-almost surely, it could be that the Gibbs measure has 
most of its mass concentrated on the ko spin configuration with ko that does 
not depend on N. We will show that this is not the case in general. Results 
for the approximated Gibbs measure and the Gibbs measure are given in 
Section 3. 

Another question is related to the /^-dependence of the Gibbs measure. 
As mentioned above, the system is frozen at the level of the free energy. 
Does a similar fact hold for the Gibbs measure? When the limit /3 f oo is 
taken after the limit N j oo, the Gibbs measure converges to a point mass 
at the spin configuration that realizes the absolute minimum as is proved in 
Section 3. So to exclude the possibility that the system is frozen at the level 
of the Gibbs measure, we estimate from below the total variation distance 
between the Gibbs measure at finite (3 > (3 C and its limit when (3 f oo. 

Assuming that some or all of the foregoing questions have been resolved — 
in particular, that the number k(N) has been determined to obtain, P-almost 
surely a good approximation of the Gibbs measure — then is there an "easy" 
way to construct this approximated Gibbs measure? This is the subject of 
the last section. 

2. Some results involving order statistics and sums of spacings. 

2.1. Notation and recollections. In this section, (X{,i E {l,...,n}) are 
i.i.d. standard Gaussian random variables with distribution function $, 
(Ui,i € {1, . . . ,n}) are i.i.d. U(0, 1), random variables and (Wi,i G {1, . . . 
are i.i.d. exponentially distributed r.v. £(!)■ We define So := and S m := 
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m G N*. In general we denote by Yi >n < ■ ■ • < Y n ^ n the order statistics as- 
sociated to some random variables (Yi,i € {1, . . . , n}). We set C/o,n := and 
U n +i,n := 1- 

Let G be denned by 1 — $(G(u)) = it, < u < 1. Then G satisfies (see 
[14], page 818), as it J 0, 



(2.1) G(u) 



;( ^ 



loglog(l/it) + log(47r) 



+ o 



(loglog(l/u)) s 

(iog(i/u)) 3 / 2 y 



2V21og(l/«) 

As it is standard, we can construct the Gaussian random variables by using 
Xj = G(Ui), and since G is decreasing, we have X^ n = G(f7 n _j + i in ). Then, 

by symmetry, we have the identity in distribution (denoted =) 

(2.2) X ijn -X ltn ^G(U 1>n )-G(U ijn ). 
Recall the results 



(2.3) 
and 

(2.4) 



{U itn , < i < n + 1} = {Si/S n+1 ,0 < i < n + 1} 



Ui. 



; 1 < i < n 



are i.i.d. U(0, 1) random variables, 



which can be found in [15], Lemmas 1 and 2, or in [43] and [32], respectively. 

2.2. New results. Before stating and proving the main result of this sec- 
tion, Proposition 2.2, we present some preliminary results: 

Lemma 2.1. For all 5 > 0, there exists hq = no(5) such that for all n > 
no, we have 



(2.5) 



log 



1 



logn 



< 2 log (logn 



,1+5 



> 1 



(logn 



,1+<T 



Proof. Equation (2.3) implies that 
1 



log 



logn 



< 2 log (logn 



(2.6) 



>P[|log5i| <log(logn) 1+5 

Sn+l 



log 



n+l 



+ log 1 + 



,1+5 



1 



n 



> log (logn 



,l+<5 



Since Si is exponentially distributed, we obtain 

(2.7) P[|logiSi| < log(logn) 1+5 ] > 1 - exp{ -(logn) 1+5 } 



(logn 



,1+<T 
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By using the exponential Markov inequality and the fact that (1 — x) 1 < 



(2.8 



for \x\ < 1/2, we get that VO < e < 1, Vn > 1, 



< 2e- n£2/ \ 











>e 




n 





which implies after a short computation that, V0<e<l, V<5> 0, 3n £)< 5 
exp(2/(l - e)) 1 /^ 4 " 5 ) such that Vn > n EjS , 

S n +1 



(2.9) P 



log 



n+1 



1 + - 

n 



> log (log n 



,1+5 



< 2e" ne / 4 . 



Then combining (2.6), (2.7) and (2.9) and taking e = 1/2 entails the exis- 
tence of no = ng(<5) such that (2.5) holds. □ 



Proposition 2.1. We have 
log 



(2.10) 



i=l 



00 'Wi 



i=3 



E[Wj} \ ^ E[Wj] 



8=1 



where Z is a random variable such that \/ x > 0, 

(2.11) P[Z > X] < e -(^/2)+15/4 mrf p[ Z <_ x ]< e -(-v / 3/2)+15/4 > 

Moreover, for all j > e, we have 

exp(3((logj) 2 /j)) 



(2.12) 



log 



Ul,n 



' 3 log j 



logj 



> 1 



,■3/4 



More generally, Vc > 0, VO < e < 1, 3 jo = jo(c,e) and V j > jo, 



(2.13) 



log 



Ujn 



<-Z + 



'clogj 



logj 



> 1 



1 



;(l-e)c/4 ' 



The proof of this proposition necessitates the following result: 

Lemma 2.2. For all positive x, for all positive integers jo, ji (with j\ > 
jo) and for all positive t such that t 2 / j" 2 < 3/4, we have 



(2.14) 



a m - E[w t ] 

> > x 

, 1 

1=30 

t 2 



< exp<^ -tx + — 1 + 2^2 1 + — 



Jo 



Jo 



jo 
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In particular, for t 2 < 3/4, 
(2.15) 



1=1 



I 



< ~tx+2t 2 (l+2t 2 ) 



PROOF. The proof is a consequence of exponential Markov inequality 
and Jensen inequality. Indeed, we can write 

* Wr- E[Wi] 
1^ 1 > x 

1=30 

<e~ tx ft E(exp\t 



W l -E[Wj] 



<e 



-t.v 



ft ^(exp{ 

1=30 V L 



=30 

*- tx ft r 

1=30 



W t -W{ 

I 



t 2 /F 



2t 

< ex.p{ —tx + I H 7T 

Jo 



1=30 > 



with (W/) i.i.d. standard exponentials, that are independent of (Wi) and by 
using the inequality (1 — x)" 1 < e x+2x2 for < x < 3/4 at the last step. Now 
(2.14) follows since 



.a 



1=30 



1 




P'i dx 






p 








v 1 + io) 


Jo 


ho x 2 







□ 



Proof of Proposition 2.1. We have, by using (2.4), 



log 



(2.16) 



Uj t n 



log 

3-1 



U 



l.n 



U 



3~l,n 



U2,n Uj t n 



1=1 



i+l,n 



3-1 



8=1 



from which we deduce the first equality of (2.10). Now we can write 



3-1 



i=l 



i=l 



Wj - E[Wj] y, Wj- E[Wj] 

i f-f i 

i=j 



3-1 

+ E 

i=l 



E[Wi] 
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where the two series converge not only in quadratic mean, but also almost 
surely by the Levy theorem. 

Calling Z = E£i(Wj - E [Wi])/i, the tail of the distribution of Z de- 
creases exponentially, as mentioned in (2.11). Indeed, it suffices, on one 
hand, to apply (2.15) for t = v3/2 to get the first property on Z and, on 
the other hand, to note that since 



i i=i 



i 



li=i 



I 



by making exactly the same computation as in Lemma 2.2, we get the same 
estimate and therefore the second property on Z. 

By applying (2.14) with x = \J{c\og j)/j and t = y/cj log j/2, we can say 
that, for all c> 0, for all j > 1 such that 3j/ logj > c, 



j^ Wj-Em > j clog j 



1=3 



I 



3 



<exp<{ --\ogj\ 1 



1 + clogj 
3 



Therefore, for all c> 0, for all < e < 1 and 3 jo = jo(c, e) such that j > jo, 
we have (1 + clogj)/j < e, which entails (2.13), whereas for j > e, taking 
c = 3 gives (2.12). □ 

Now we can state the main result of this section. 



Proposition 2.2. For all < 5, for all < e < 1 and for all k n satis- 
fying 



kr, , logk n 
j oo and I (J 



(2.17) A: n Too, 

log n log n 

by defining, for < A < 1 and < a < 1, 

logn 



as oo, 



(2.18) 
and 
(2.19) 



An, :— A- 



X n := X n 1 + 2 



2(1 -A ra )' 

AnP-")/ 2 



there exists no = no(e, S, A, a) smc/i i/iai Vn > no, i/iere exists fi n C 0, wii/i 



(2.20) 



[fij > 1 - 4 



1 



-kn/16 



+ 



(l g n ) 1+5 !_ e -l/16 

such that on Q n , for all j such that k n <j < nX n , we have 



(2.21) V2 log n{X hn - X 1>n ) > 2 log n 1 - J 1 



logn 
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while for all j such that nX n < j <n, 

(2.22) V21ogn(Xj, n - X hn ) > 21ogn^/T^e - G(A n )V21ogn. 
In particular, on £l n , for log j / log n f 1, 

(2.23) V21ogn(X jjn - X hn ) > 21ogWl - e - G(A„)\/21ogra. 
Remarks. 

1. If A n is chosen as a constant A independent of n, < A < 1, (2.22) gives 
that Vj > nX, y/2logri(Xj tn - X 1)Tl ) > 2\ogny/l -e - G(A)( v / 21ogn) > 
21ogn(l — e) if n > no(X,e) for some no(A,e). When entering various 
regimes as nA^ < j < n\n with An I for i = 1, 2, a cancellation could 
occur between the two terms in the right-hand side of (2.22). The choice 
of X n in (2.18) then allows us to see such cancellation, since for a < 1/4 
it provides that, Vj > nX n , on Q n , 

(2.24) v / 2logn(X J> - X hn ) > 2(1 - y/a)(l - s) logra > (1 - e) logn. 

2. Note that the lower bounds in (2.21) and (2.22), even when obtained by 
two completely different methods, are of the same order 2(1 — \fa ) logn 
when j = nX n . 

3. Note that for n = 2 N , under (2.17), we have J2n=i f [^ c 2 n] < oo. In the 
next section, this allows us to get some results that are true P-almost 
surely for all but a finite number of indices N by the first Borel-Cantelli 
lemma. 

4. For completeness, even though it will not be used in the next section, we 

show that for < rj < 1, there exists Q* kn with P(fijj^) > 1- 2 exp{- V3/2(log k n ) v }, 
such that on Q% , for all j such that k n < j < nX n and log j'/ logn [ 0, 
(2.21) can be refined as 



Proof of Proposition 2.2. Let us separately consider the two cases 
as in (2.21) and (2.22). 

Case (i). j is such that nX n < j <n. Let us denote by 



(2.25) 




n 



(2.26) £ = £(n,A n ,e):= (J {v / 2bg^(X n , n - X, 



j+l,n) < f(n,X)} 



j=n\„ 



with f(n, A) := 21ogn-y/l — i — y / 21ognG(A n ). 
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Note that in this case we are working where j is large with {X n ^ n — 
X n _j + i :n ) instead of (Xj iU — X± in ) since these spacings are equal in distri- 
bution. We claim that 

(2.27) P[£] < 2e~^ 1 ~ x ^ n + F[X n<n < y/2(l - e) logn], 

which entails that VO < A < 1, 3 no and Vn > no, 



(2.28) 



< -t + F[X n , n < V2(l-e)logn], 



71- 



after having noticed that VA, < A < 1, Va, < a < 1, 3 no = no(A,a), s.t. 
Vn>n , 4(1- A„) > 2. 

Now, (2.27) is an immediate consequence of 

(2.29) F[C n {X n:n > V2(l-e)logn}] < 2 e - 4 ( 1 - A ") lo s™. 

To prove (2.29), note that, on one hand, 

(23Q) £n{X n , n > V2(l-£)logra} 



C£:= |J {V21ogn(V(l - e)21ogn - X n -j+l,») < /(n, A)} 
j=nA„ 

and, on the other hand, 



(2.31) F[C n (J5f n , n > V2(l -e)logn )] < F[A] + F[£ n -4 C ] 
with 



4:= 



(2.32) 



i=l 



X^-p^G^n)} ~ E [ 1 {X l <G(A n )}]) 

>A n ,(l-A n )n (1+a)/2 2V27Al, 



where E[l{x. i <G*(A n )}] = 1 ~ by definition of G. 

To estimate P[»4], we use Bernstein's inequality [41], namely: 

Lemma 2.3. Let {Yi)% be independent random variables, mean 0, s.t. 
E[Y?} <oo and \Yi\<l,Vi. If 



D n :=J2E[Y?] 



i=i 



and if < t < y/D n , then 



i=l 



>WD n 



<2e^ 2 / 2 . 
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By usin g this inequali ty for Y { := t{ Xl <G(\ n )} ~ ^{Xi<G(X n )}] with 
t = 2n a / 2 y/\ n (l - A n )2/A, since D n = nA„(l - A„), 3n (A,a), Vn > no, 
t < \/D n , therefore we get that V n > uq , 

(2.33) P[A] <2e" 4 ( 1 - A ") logn . 

On the other hand, 



(2.34) A c C n(l - A n ) - A n (l - X n )n^/ 2 2V2/X < £ t {Xi 



{Xi<G(X n )} (, 



i=l 



which implies that on A c the number of random variables Xi less than G(X n ) 
is greater than 1 + nil — A n ); hence, A c C {^i +n (i-\ n ) n < G(A n )} and so 



(2.35) 



Sn—nXn+l \ 

(J A\ n > G(A n ) n (A\ 
^ fc=i / 



+n— nX n ,n 



<G(A n )) 



0. 



Combining (2.31), (2.33) and (2.35) entails (2.27). 
Now (2.28) and the fact that for all < e < 1, 

(2.36) F[X n>n < v / 2(l-e)logn] = $ n (V2(l-e) logra) < e -«7V 2 ( 1 -<0 1 °s fl 
give 



(2.37) 



n 2 



which leads to (2.22), since, for all < e < 1, for all 5 > 0, for all < A < 1, 
for all < a < 1, for all /c n satisfying (2.17), 3 no = no(e,5,\,a) such that 
Vn > no, 



(2.38) A + e -^/V 2 ( 1 - £ ) 1 °g«<2 
n 2 



1 



-fc„/16 



+ 



(l ogn )i+5 i_ e -i/i6 



Case (ii). k n < j < n\ n . Note that, because of (2.1), we have, for u J, 0, 



G(u)-G(u) = V2^1ogQ 



'log! i 



+ 



log log (1 /it) 



loglog(l/t>) 



For 7 > 0, let 



n X» == { 



sup_ U j>n = U n - Xntn < (l + 7)Ar 
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Since 

mi ] < 

An 

< 
< 



n\ n 

An 



S 



n+l 



n + l 



1 



< e 



+ P 



S. 



11+1 



n + l 



1 



S 



n\ Tl 



Sn+1 

s. 



>(l + 7 )A n „5 n+ i>(n + l)(l-e) 

"n+l 



nA ^>(l + 7)(l-e) n 



nX 



n 



n + l 



> £ 
+ P 
- 1 



n+l 



n+l 



>£ 



by using (2.3) in the second inequality, it follows from (2.8) with e = 7/(2 + 
I } < 4exp(-nA n ( 
ie have 

G(U 1>n ) - G(Uj, n ) 



7), that F[Ql ] <4exp(-nA n (7 2 )/(4(2 + 7 ) 2 )) 
On .we have 

An 



(2.39) 



+ 

Let Q n < 5 C 0, be defined by 
Then Lemma 2.1 gives 

(2.40) p[n n , s ] > 1 



loglog(l/A n )\ 



Vlog(l/A n ) 



logn 



< 2 



log (logn 
logn 



,1+<J 



(logn) 1+<5 ' 



Let r2 n>< 5 := ^l U) s H . Then by combining (2.2), (2.5) and (2.39), on Q njl 
we obtain that 

V2logn(Xj tn - X 1 . n ) 

(2.41) 

= 2 log n 



Again using (2.3), we can write 

log(^i, n /C/ 3> ) _ log(Si) logQVi) log j 



(2.42) 
Let 

(2.43) &j :-- 



log n 



log n log n log n 



< - [ and ft kn := I f| %j fl O n ,5- 
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Then by using (2.40) and (2.8) with e = 1/2, we get 



(2.44) P[^J<4 
On we have 
(2.45) 



1 



-fen/16 



+ 



(l ogn )l+<5 1 _ e -l/16j- 



log(Ul,n/Uj,n) _ log .7 



log n log n 

which combined with (2.41) gives 



+ o 



log ( log n 
logra 



,1+5 



(2.46) 



21ogn 1 - Wl 



logi ^ log (log n)^ 



logra 



logra 



\ Vlog?7, 



Note that by Taylor's expansion, the right-hand side of (2.46) is of the form 



logjl l + O 
Hence, if we define a n j by 



log log n 
log fe n 



l + O 



log log n 
ydogn 



a n j log j := (l-e n )21ogn 1 - Jl 



logj 
log re 



with e n well chosen, then we get 



I (J V2logn(Xj >n - X 1>n ) < a n j log j J n Q k 

-\k n <j<n\~ 



and 



(2.47) 



f (J V21ogn(X 3> - Xi, n ) < a nJ log j J n VL% n 



'fen<i<nA n 

<P[fi£J<4 
and therefore (2.21). 



-fc„/16 



+ 



(l ogn )i+<5 i-e-i/iey 



To prove (2.25), that is, when j 6 Jfc n := {j £ [fe n ,nA n ], such that 
log j/logn | 0}, by Taylor's expansion in (2.41), we have on Cl nj s, 



(2.48) 



v / 2logn(X 3> - X 1>n ) 



( Ul,n 
\Uj,n 



log 



+ o lo 



Ul,n 



l + O 



log log n 
V / Iog~n 
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Therefore, to evaluate 



(J v / 21ogn(X J> - X hn ) < a* „ log j n n nj6 



where a* := 1 — l/(logj) 1 ^5 it is enough to estimate 



k 



U -loi 



.( Ul,n 

'\U j>n 



< a* jiV iogj nn n ^ 



By using (2.13) with c such that (1 — e)c/4 > 1 and jo = k n , and by defining 



Ml 



n Jog 

j>k n K 



U; 



< -Z + 



we get that P[fi£j > 1 - 

p fc „<p U \ z < ( 



(l-e)c/4-l. 



/clogj 



logj 



. Now it is easy to check that 



jlogj 



+ (a* >r) — 1) ) log j 



<¥[Z<-\x kn \logk n }+n(n* kn ) c } 



with \xj \ := l/(logj) 1 v — \J cj j log j, where we have used at the last step 
that (\xj\ log j)j>k n is an increasing function of j. Hence via (2.11), we get 



V3, 



(2.49) P fcn <e X p|-^-(log/c n ) r ' + ^- 
which implies (2.25). □ 



\/3 c\ogk n 15 



+ ^ + 



(l_ e )c/4-l 



Proposition 2.3. For all S > 0, /or a// fc n satisfying (2.17) and /or 
C giwen fry fi* := 0^ n O nj 5 n Q kn defined, respectively, in (2.52), 
(2.40) and (2.43), we /iaue ' 

(2.50) p[o;] > 1 



and on Q n , we have, V j, 1 < j < k n , 



(logn 



Wi 



,1+5 ' 



(2.51) v / 2bg^(X i!n -Xi, n ) = ^^- + 



i=l 



log (kn) log logra 
logn 



Proof. For 7 > let us define 
(2.52) n n>k :=\ sup U j>n = U kn>n < (l + 7 ) — )• 
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It follows from (2.3) and (2.9) that P[fi£ >fc ] < 4exp(-/c n (7 2 /(4(2 + 7 ) 2 ))). 
From now on, our work space will be & n: k- By using (2.1), we can write 



V21ogn(X, in - X 1>n ) 



a/2 logra. / 2 log ' 



L / \Qg{Ul,n/U h n) 

I V io g (iM, n ) 

loglog(l/^i, n ) | loglogCl/^-, 



4Iog(l/^i, n ) 4 v /log(l/f/ J , ri )^log(l/C/ 1 , n ) 
log(4vr) / 1 1 



^(^l,.) Vlog(l/C/ J> )^/log(l/^i, n ), 

- loglQg(lM, n ) \ 2 | (l0gl0g(l/L/ Jin ,)) 2 



io g (i/cr lin ) ; iog i / 2 (i/ C/lin )iog 3 / 2 (i/c/ J - n );j' 

but by using the same type of arguments as in the proof of Proposition 2.2, 
Case (ii), we have, on f2 n ^nOfc n [defined in (2.40) and (2.43), resp.], 

1 1 



sup 



l<i<fcnVl0g(l/C/l, n ) Vl0g(l/^, n ,)^l0g(l/C/ lin ) 

1 



log(l/C/ 1)n ) V v /(log(l/C/ fcn , n ))/(log(l/C/ 1 , n )) 
log k. 



O 



log 2 n) ' 



sup 



loglog(l/C/i, n ) loglog(l/C7i in ) 



V log(lM,„) Vlog(l/[/ j , n ) v /log(l/C/ 1 , n ). 
_ loglog(l/[/ 1 , n ) / 1 \ 

iog(i/J7i n ) V .Awn iu,. _^/rwn 7r^ 



= o 

and also 



log(l/E/i, n ) V v /(log(l/C/ fcn , n ))/(log(l/C/ 1 , n ))' 
log(fc n )loglogn N 



log 2 n 



doglog(l/[/ 1 , w ) \ 2 | (loglog(l/[/ iin )) i 



i<3<k\\ log(l/^i,„) / log 1 / 2 (l/[/ lin )log 3 / 2 (l/C/ iin ) 
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log(/c n )loglogra s 



log 2 n 

Therefore, we obtain that, on Q n> g n Ofc n , Vj, 1 <j < k n , 



y/2]ogn(Xj n - Xi, 



o /, , 1 J, /-, . l°5(Ui,n/U jtn ) , /log(/c n )loglogn 

2 J log n log—- ^ 1-W1 + — — ' ■" +0 — g 

V U ln I V log(l/i7i,„) V log 2 n 



/log n log J- [ - MgW^) + ( log(fcn)loglogn 



C/i, n l log(l/E/i, n ) V log 2 n 



^2~ + 0^ 



i=l 



z V log n 



by using Taylor's expansion in the second equality, and (2.10) and (2.40) in 
the last equality. □ 



3. Applications to REM. We are now able to answer the questions raised 
in the Introduction by using the tools developed in the previous section. We 
choose n := 2 N and use N instead of n(N), since N is related to the volume 
of our system. Otherwise, we keep the same notation used up to now, in 
particular C = a/2 log 2, with 0/(3 c > 1 . 

The Gibbs measure, defined in (1.4), can also be expressed as 



2 V 2 tfffffVIe-^* 

(3.i) w (*) = £ n(*)i)mA(°)i) = l % 2 , _ pVNXi — . 

i=l 2^i=l e 1 

where (X { ,i £ {1, . . . ,2 N }) are i.i.d. standard Gaussian random variables [as- 
sociated to a = {(a)i,i E {1, . . . , 2^}}, an enumeration of { — 1, +1} ] and 
is a real function defined on { — 1, +1}^, such that H^Hoo := sup CTg {_ 1)+1 }jv |^(o")| < 
oo. 

We denote by (a)j the (a)i to which Xi = Xj 2 N is associated, so we can 
write 

„ „.,, E?.i1'(W<)^P(-Wft)v / 21og2' v X,^) 
( ' W< ' = Eg,ex P (- Wft )V2 1 o g 2^„ 2 „) • 

3.1. Approximation of the Gibbs measure. The first result of this section 
is related to how many terms of the order statistics we need to ensure a 
good approximation of the Gibbs measure. 
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Theorem 3.1. Let satisfy 

(3.3) kjy | oo, 1 oo and — ^ — [ asN^co. 

There exist Qjy C SI and /Vo siic/i i/iai for all N > No, 

and on Oat we /iaue 



(3.5) _ ^((erji) exp(-(/3//3 c )v / 21og2 Ar (A^ 2J v - X^ N )) + 



Ei^M-(P/f3cW2log2 N (X it2N - X lt2N )) + B N (1) 



with 



(3.6) : 



P/P c -l{k N -\)WP°-W 

PROOF. From (3.2), we deduce (3.5) with B N (V) := Ei=k N +i x 
exp(— (0/0 c )V2 log 2^ (X i2 N — X 12 n)). We have to prove only that sat- 
isfies (3.6), which is a consequence of Proposition 2.2. 

Let ifN(i) '■= exp(— (/3//? c )v / 21og2 Ar (A" ij2 A' — -^1.2^)) an d := ^2 N 

satisfy (2.19). Then 

/2 N \ N 2 N \ 

Bjv(tf)<||*||oo( S ¥>i\r(«)+ X] 
On one hand, by using (2.24), we have on Qn, 



2 N 

£ MO < expj- (£(1 - e) - l) k^j = T 



Nt 



i=2 N \ N +l 

(3.7) 



with r := -£-(1 - e) - 1 >0. 
On the other hand, by using (2.21) with e = (1 — c /0)/2, we have, on Jl n , 
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E m(i) 

k N <i<2 N \ N 



E exp(-|-(l- £ ) , 2lQg ^ } 

(3.8) 



k N <i<2 N \ N 



< 



E 



tCl+rXl-e) 
k N <i<2 N X N 



Hence (3.7) and (3.8) imply (3.6). □ 

Remark. An example of k^ that satisfies (3.3) is k^ = Nlog p (N), 
where log p = loglog p _ 1 . 

Now we can define the random probability measure on {— 1, +1}^ by 



(39) (i) m . E^i^((^ fc )exp(-(/?//3 e )V21og2^(^ 2J v -X li2N )) 

and thus show that this measure is a good approximation of the Gibbs 
measure. Indeed, by using the total variation distance between two measures 
given by 

(3.10) <hv(ji,v)= sup 

{*: ||*||oo=l} 

the first Borel-Cantelli lemma and the fact that the denominator in (3.5) 
can be bounded from below by 1, it is immediate to check the following 
corollary: 

Corollary 3.1. For all (3 such that (3/[3 c > 1, for all k^ that satisfy 
(3.3), we have with probability 1, for all but a finite number of indices N, 

(3.11) dMm,p^i l lp) < p/ l-i {kN -i)W^-m- 

With this last result and Proposition 2.3 in hand, we can define a second 
approximation of the Gibbs measure by the random measure on {— 1, +1}^, 

(2) f\vt\ ._ E k k ^mk)eM-P/PcZe=iW e /£) 



(3-12) tfijm : 



Eti^M-P/PcELiWe/i) 



where (Wi, 1 < I < kjy) is a family of i.i.d. exponential random variables 
with mean 1. 
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Proposition 3.1. There exists an absolute constant C such that, for 
all f3 satisfying f3/(3 c > 1, for all kjy satisfying 

/q i n\ / * kjv + , (logfcjy)(logJV) 

(3.13; kn\oo, — \ oo ana — I (J aslV^oo, 

we have with probability 1, for all but a finite number of indices N: 

dTv{fi<k N ,f3,lJ<kl t p) 

,« u n < r { (logk N )QogN) \ f (logA^XlogiV) ) 

( J " A I N J eXP l C /3 c iV / 

4 1 

+ 



P/Pc-likN-qWPc-W 
vc 

rf Tv(^ l/3 ,/4 2 2,/3) 



Proof. We have 

"TV V> 

(3 ' 15) <c /3 / (log (log jV) \^ p f c /? (logfc JV )(logJV) 



Pc\ N J { f3 c N 

as a consequence of Proposition 2.3, the fact that Wx, \e x — 1| < |x|e' x ' and 
the first BorePCantelli lemma. Applying the triangle inequality, (3.11) and 
(3.15) gives (3.14). □ 

Remark. Note that (3.13) is satisfied if we choose k^ = N\og p (N), 
where log p (x) = loglog p _ 1 (x) and log 1 (x) =logx. 

3.2. Some properties of the Gibbs measure. The previous results help us 
to get more information on the support of the Gibbs measure. Indeed: 

Property 1. As an immediate consequence of Corollary 3.1 and Propo- 
sition 3.1, we have: for kjy satisfying (3.13), 



(3.16) limsup/iAr i/3 ({(a)i,...,(cr) fcjv } c ) = a.s. 

JVtoo 

Property 2. When /3 C < (3 < oo, the Gibbs measure cannot be concen- 
trated on the minimum, that is, ^k N ,p 7^ ^(o-)i- 

Proof. Considering (3.12) with the trial function ^((a)i) = lr i=1 tVi > 
1, we can write 

d Tv(4 2 j,/3>%)i) 



GIBBS MEASURE FOR REM 21 
> Efeexp(-/3//3 c Eii^/^) 



> 



l + Efe 2 exp(-^//? c EiiW,/£) 
expC-jS/^Wi+Wz^)) 



l + exp(-P/0 e (Wi+W 2 /2)y 
which in combination with Proposition 3.1, provides that 

- l + exp(-^ c (W 1 + W 2 /2)) >0 - □ 
To go on with the properties, we need the following lemma. 
Lemma 3.1. For all j3 > 0, 



a.s. 



(3.17) lira £exp £Wi//) = C(/?) 

If f3> /3 C , then E [((P)] <oo. In particular, P[0 < (((3) < oo] = 1. Moreover, 
for any given sequence (3 P such that limp-foe /3p = oo, i/tere exists a subset 
Q C f2 snc/i i/iai P[il] = 1 and, on fi, 

(3.18) hmC(/? P ) = 0. 

pjoo 

Proof. Let £fc be the a-algebra generated by (Wi, . . . , Wk) and let 

The positive random variable Ck(P) is a super martingale, since we have 

(3.19) E[C fc+ iC9)|E fc ] = — — J 7 — — CfcC9) < C*C9). 

1 + + 1 )) 

Therefore, the increasing sequence Cfc(/3) converges almost surely to some 
random variable C(P)- On the other hand, using that Vx > 0, log(l + x) > 
x - (x 2 )/2, ELiW" 1 > lo §i and ELiW" 2 < 2 - 1/i < 2 provides 
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where at the last step we have used (3/f3 c > 1. Hence, by the Fatou-Lebesgue 
lemma, we get E [£(/?)] < oo and so P[0 < ((/3) < oo] = 1. 

Proving (3.18) requires a little more care, since the probability subset 
where Ck(P) converges, denoted by £l(f3), depends a priori on (3. 

We can write ( k (/3) = exp(-(J3/(3 c )Wi)( k {P) with 

cm : 

By the same kind of computations we did previously, we can show that, 
for any sequence (3 C < (3 p f oo, there exist some positive integrable random 
variables (((3 P ) such that on the subspace Cl := f)p fi(/3 p ) of P probability 1, 
we have lim fcToo ( k {/3 p ) = C{P P )- 

Since for all k, Cfc(Aj) < Ck(P P ) when (5 p < q , we have, also on Cl, Ciflq) — 
C(j3p). Therefore, denoting f3* = inf{/3 p > l + p e }, since P[C(/?*) < oo] = 1 and 
Pflim^-foo exp(— ((3/(3 c )Wi) = 0] = 1, with probability 1 we get that 




< limsu P C(/3 P ) < limsupe (_/3//3c)Wl C(/3*) = 0. 



Property 3. At zero temperature, the Gibbs measure is a.s. the point 
mass at the minimum. 

PROOF. We have 

1 + E fc =2 exp(-/3//3 c ELi 
which gives with probability 1, via Lemma 3.1, 

TVtoo 



limsupd T v(/4v,/?>^)i)^ l + apY 



Then Proposition 3.1 provides with probability 1 that 

2C(/3) 



limsupd T v(wv,/3,%r)i< 



< 



Now by using (3.18), we get with probability 1, 

limsuplimsupdTv(A i Ar,/3) <^(ct)i) =0. 



□ 



Remark. Note that there is a result in [29] , in the case of "REM for size 
dependence," that deals with the quality of the concentration of the Gibbs 
measure on the minimizer of the energy when the temperature goes to zero 
with the system size. 
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Property 4. With probability 1, the Gibbs measure is not the uniform 
measure on the k^ first minima {(cr)i, . . ., (&)k N }- 

Proof. Let v^ N be the uniform measure on kjy points. Considering 
(3.12) with ^((cr)j) = l(j = i) Vi > 1, we can write 

> 1+Et!eip( ^ ELw() - 

which, combined with Lemma 3.1 and Proposition 3.1, gives that with prob- 
ability 1, 

liminf d TV ( W; y kN ) > > 0. n 

Property 5. With probability 1, the Gibbs measure is not the measure 
Mfc ,/3 for a finite k . 

Proof. Because of Proposition 3.1, it is enough to prove that for a 

(2) (2) 

finite ko, lim inf jvtoo rfTv(/-4 a-, \^\ N r) is bounded from below by a quan- 
tity that does not go to zero. Considering (3.12) with ^ defined by ^{a G 
{(<5")i, • • • , (o")fc }) = an d = = 1 for all ko < i < ku, we get 

f3 21) d ^^klp^kl,p) > m£,/?({(^)i, • • • . (°) ko Y) 

~ i + E^ 2 exp(-m£f=iW^)' 

and so, via Lemma 3.1, with probability 1, 

v • fr , / (2) (2) w exp(-^ e E^^) . n 

Property 6. With probability 1, for j3 c < (3 < oo, the Gibbs measure 
does not have all its mass concentrated on {(cr)i, • • • , (cr)fc } f or a finite k$. 

Proof. This result comes from ^k ,/3[{(o')i, ■ ■ ■ , (cr)fc }] = 1 an d 

//A^[{(cr)i,...,(a-)fc }] 

= MJV,/3[{(c0l, • • • , - Mfcjv^K^)!' • • • > (^)fcol] 

+ Pk N ,p[{(*)li ■ ' ■ ' (^feol] - ^fco,/s[{(ff)l, • • • , (^)fco}] + 1 
if one notes that 

MfcAr,/?[{(^)i> ■ • ■ > (^)feol] - Wo,/)[{(^i. • • • > (&)ko}] 

= Mfco,/3[{(ff)i, • ■ • , (<7)fc } C ] - ■ • ■ > Fl 

= -Mfe w ,/3[{(^)l»--->(^)ifco} C ]- 
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By using (3.21) and Proposition 3.1, we get 

limsup/ij Vj a[{(5-)i, . . . , (cr) feo }] 

Nloo 



k +l 



3.3. Representation of the Gibbs measure. Thanks to the two approxi- 
mations of the Gibbs measure given in (3.9) and (3.12), we are now able 
to propose a representation of the Gibbs measure and a way to simulate it. 
First we consider the measure /ijjj^ p defined in (3.9), that we can simulate 
in the following way. [Note that we have to pay attention to the relation- 
ship between the (<7j) and the (uj), since the (<7j) are needed for the simula- 
tion.] Let (Ui, . . . , U 2 n) denote a 2 N sample of uniformly distributed random 
variables of order statitistics (Ui 2 n < • • • < U 2 n 2 n). We want to construct 

exp{— ((3 / (3 c )y/2\og2 N (X k 2 n — ^1,2^)}' because of (2.2). We consider and 
simulate only the k^ last terms U 2 n 2 n, . . . , U 2 N_ kN+l 2 n, with /cat satisfy- 
ing (3.3). (Note that this can be done without having first to order uniform 
variates by using some faster algorithm (cf., e.g., [23], Section 3.26).) Then 
independently, we choose, one after the other, spin configurations in 
{— 1,+1}^ without replacement. This defines an ordered sequence of con- 
figurations that we call ((c)i, . . . , (a)k N )- We make the following claim. 

Claim 1. The random measure defined by 




(3.22) 

+ £ vi;(( (J ) fc )exp(-(/3//3 c )v / 21og2 Ar (G(f/ li 2iv) - G(U k>2 „)) 

k=2 ) 
( k N \ -1 

x 1 + E exp(-(/3//3 c )v / 21og2 7V (G(C/ li2 iv) - G{U h ^))) 

I k=2 ) 

has the same distribution as a. 

Indeed, the claim is an immediate consequence of the following technical 
lemma on spacings of n-independent uniform random variables. For 1 < j < 
n, let £j be the index of the (almost surely unique) Ui such that Ui j = Uj >n . 
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Lemma 3.2. For all (xi)i<i<k and for all integer k, we have 
W[U i:n <x i ,U i , n = U li ,l<i<k] 

= W[U ijn <Xi,l<i<k] x P[U itn = U h ,l<i<k]. 

Moreover, 

in - kV 

PROOF. We have (cf., e.g., [13]) 

W[U ijn <Xi,l<i<k) 

u 
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r 3j i p x k 

/ ••• / ( l -yk) n ~ k Ho<yi<y2<-<y k <i)I[^- 

i=i 



(n-k)\ 

Moreover, we can show by induction that 

W[U iin < Xi ,U i:n = U k ,l<i<k] 

rxj rx k k 

= ■■■ (1 " yk) n ~%0< yi <y 2 <-<y k <l) II d Vi> 
JO JO i=1 

which completes the proof. □ 

Thus we have the first representation of the Gibbs measure that can 

(2) 

be simulated. An alternative and easy way to proceed is to consider fi k ^ a 

defined in (3.12) instead of //jj. p. To simulate ^ a, we just have to con- 
sider a kN sample of independent uniformly distributed random variables 
U\, . . . , Uk N , with k n satisfying (3.13), then choose as before the kjy spin 
configurations, (cr)i, . . . , (cr)k N to construct the resulting measure 

(3 23) u (2) m - ^^ + ^=2^^)eM+(3/(3cEi=i(^Ui)/^ 
kN ' P l + E^ 2 exp(+/?//3 c Ef=i(log^)/^) 

(2) rl (2) 

Then we have, in the same way as for the above claim, that jl k ^ a = yr k ^ p. 

Thus this second procedure needs only two independent samples: one of 
kpf spin configurations chosen without replacement within 2 N and one of fc/v 
independent uniform random variables on (0, 1), so we do not have to work 
with the k n largest order statistics. 

Note that, as mentioned previously, k^ must satisfy (3.13) and so could 
be chosen as N log log log N, for instance. 
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